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ABSTRACT 

Aims. We present a study of the evolved stellar populations in the dwarf spheroidal galaxy Fornax based on wide-area 
near-infrared observations, aimed at obtaining new independent estimates of its distance and metallicity distribution. Assessing 
the reliability of near-infrared methods is most important in view of future space- and ground-based deep near-infrared imaging 
of resolved stellar systems. 

Methods. We have obtained JHK imaging photometry of the stellar populations in Fornax. The observations cover an 
18.5 x 18.5 arcmin 2 central area with a mosaic of SOFI images at the ESO NTT. Our data sample all the red giant branch 
(RGB) for the whole area. Deeeper observations reaching the red clump of helium-burning stars have also been obtained for a 
4.5 x 4.5 arcmin 2 region. 

Results. Near-infrared photometry led to measurements of the distance to Fornax based on the K-band location of the RGB 
tip and the red clump. Once corrected for the mean age of the stellar populations in the galaxy, the derived distance modulus 
is (m — M)o = 20.74 ± 0.11, corresponding to a distance of 141 Kpc, in good agreement with estimates from optical data. We 
have obtained a photometric estimate of the mean metallicity of red giant stars in Fornax from their (J — K) and (V — K) 
colors, using several methods. The effect of the age-metallicity degeneracy on the combined optical-infrared colors is shown to 
be less important than for optical or infrared colors alone. By taking age effects into account, we have derived a distribution 
function of global metallicity [M/H] from optical-infrared colors of individual stars. Our photometric Metallicity Distribution 
Function covers the range —2.0 < [M/H] < —0.6, with a main peak at [M/H] ~ —0.9 and a long tail of metal-poor stars, and 
less metal-rich stars than derived by recent spectroscopy. If metallicities from Call triplet lines are correct, this result confirms 
a scenario of enhanced metal enrichment in the last 1-4 Gyr. 
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1. Introduction 

Stellar populations in dwarf spheroidal galaxies (dSph) are 
important for our understanding of galaxy formation and 
evolution. Dwarf sphcroidals in the Local Group can be 
studied in detail, giving strong constraints on the star for- 
mation history (SFH) and chemical evolution of these sys- 
tem. Galaxies of the dSph type all started forming stars at 
an old epoch (> 10 Gyr), but in most cases this early stel- 
lar generation was later followed by major star-formation 
episodes giving rise to significant or even dominant inter- 
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mediate age populations. While old and intermediate age 
stellar populations of dSph galaxies have be en the subject 
of many st udies in optical bands (see, e.g., iGrebell 12005 : 
Held 20051 and references therein), they have been little 
studied in the near- infrared (NIR). 

However, NIR bands have several advantages when 
studying evolved stars in these stellar systems. Infrared 
photometry of evolved low-mass and intermediate-mass 
stars on the red giant branch (RGB) and the helium- 
burning phases (e.g., red clump, RC) can be used to de- 
rive the basic properties of galaxies (distance, metallicity). 
Techniques to measure such properties from near-infrared 
observations are becoming increasingly important since 
the NIR wavelength domain will be central to future in- 
strumentation (ELT adaptive optics, JWST, etc.). There 



2 



M. Gullieuszik et al.: Near-infrared observations of the Fornax dwarf galaxy. I. The red giant branch 



are advantages in using near-infrared photometry for RGB 
stars, and even more so in combining optical and NIR 
data. As we will show in this paper, the age-mctallicity 
degeneracy affecting the color of RGB stars (and there- 
fore metallicity determinations) is much less severe using 
optical-infrared colors. 

In order to explore the information contained in the 
near-infrared spectral window, we have undertaken an 
imaging study of the evolved stellar populations in Local 
Group dwarf galaxies. The first galaxy we consider is 
Fornax, one of the most interesting cases to study, be- 
ing one of the most ma ssive and lum inous dSph satel- 
lites of the Milky Way ([Mated fl998h . This galaxy was 
one of the first to provide evidence of an intermediate 
age stellar population, probed by the presence of lumi- 
nous carbon star on the a s ympt o tic giant branch (AGB ) 
(|Aaronson fe Mould! fl980l ll985t lAzzopardi et al.l Il999h . 
Another indicator of a conspicuous intermediate-age stel- 
lar population is the well populated red clump of helium - 



Table 1. Observation log. 



burning stars ([Stetson et al.lll998t ISaviane et al.ll2000al ). 

A wide plume of main-sequence turnoff stars indi- 
cates that the star formation history of Fornax has been 
continuous from ga l axy's formation up to recent times 
(IStetson et al. 19981 Buonanno et al~lll999l : ISaviane et al 



2000aUPont et alJl2004l ). The luminosity of the brightest 



blue stars shows evidence that Fornax has been forming 
stars at least up to 200 Myr ago, w hile it does not con- 
tain any stars young er than 100 Myr ([Stetson et al. I ll998l 
Saviane et al.ll200"0a ). An old population is also present, 
as shown by t he presence of an horizontal branch (HB ) 
and RR-Lyrae (|Bersier fc Woo J2002t iGreco et al.ll2006t ). 
although the blue part of the HB is poorly populated, 
so that the old and metal-poor population must be small 
( Stetson et al. 19981 : Buonanno et al.lLl999l : ISaviane et al. 
2000alt . The picture emerging from the above mentioned 



studies indicates that Fornax began forming stars at 
the epoch of formation of the Galactic globular clusters 
(GGCs) ~ 13 Gyr ago. The star formation rate was quite 
low in the first Gyr, and then increased rapidly; Fornax 
formed most of its stars 4-10 Gyrs ago. Given the presence 
of recent star formation, one would expect to find some 
gas associated with the galaxy. lYound (Il999h searched for 
neutral hydrogen out to the ti dal radius, and f o und n one. 
More recent observations by iBouchard et al. ( 2006[ ) re- 
vealed an extended Hi cloud in the direction of Fornax 
that may (or may not) be within the galaxy. 

Two studies of the chemical enrichment history of 
Fornax from the spectra of red giant stars, using the 
Ca II infrared l ines equivalent widths, have recentl y 
been prese n ted (IPont et all 120041 iBattaglia et~aTl l2006h . 



Pont et al 



(l2004r 7 found a metallicity distribution of 
Fornax centered at [Fe/ H] = —0.9 (on the scale of 
Carretta fc Gratton|[l997l . hereafter CG97), with a metal 
poor tail extending to [Fe/H] ~ —2.0 and a metal-rich 
population reaching [Fc/H] ~ —0.4. The derived age- 
mctallicity relation is well described by a chemical evo- 
lution model with a low effective yield, in which an initial 
rapid enrichment is followed by a period of slower enrich- 
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ment, reaching [Fe/H] ~ —1 about 3 Gyr ago. Then a high 
star formation rate produced an acceleration, increasing 
the [Fe/H] to a recent value of ~ —0.4 dex. These results 
agree with those pr eviously obtained for a few stars by 
(jTolstov et al. 120031) using high resolution spectroscopy. 
In a large spec t rosco pic study of RGB stars in Fornax, 
IBattaglia et al.l (2006!) noted the lack of stars more metal 
poor than [Fe/H]= —2.7 and confirmed the presence of 
a metal-rich tail up to [Fe/H] ~ —0.1. They also found 
that the metal-rich stellar populations are more centrally 
concentrated, having also a lower velocity dispersion than 
metal-poor stars. 

This paper presents the results of a near- infrared study 
of RGB and RC stars in Fornax, aimed at obtaining new 
independent estimates of its distance as well as informa- 
tion on its metallicity distribution function from a combi- 
nation of optical and NIR photometry. This study may 
represent a local example of NIR studies of more dis- 
tant resolved stellar systems with future ground-based and 
space instrumentation. Section[2] presents the observations 
and the reduction, with special care for the mosaicing 
techniques, and provides photometric catalogs. 

Color-magnitude diagrams are presented in Sect. [3] In 
Sect. U the luminosity functions of red giant stars and 
helium-burning red clump stars arc obtained and used to 
provide new estimates of the distance to Fornax. The mean 
metallicity and metallicity distribution of red giant stars 
is derived in Sect. [5] from their (V — K) color distribution 
and compared with recent spectroscopic work. Finally, our 
results are summarized and briefly discussed in Sect. [6] 



2. Observations and data reduction 

2.1. Observations 

Near-infrared observations of Fornax were carried out in 
November 10-11, 2000, using the SOFI camera at the 
ESO NTT telescope. The camera employed a 1024x1024 
pixel Hawaii HgCdTe detector which was read in Double 
Correlated mode. We used SOFI in Large Field mode, 
yielding a pixel scale of 0'.'29 pixel -1 and a total ficld-of- 
view of about 5' x 5'. 

Two complementary sets of images were taken in the 
JHK S filters. The first was a wide-area series of 16 shal- 
low contiguous fields, in a square array of 4 x 4 cov- 
ering about 18.5 x 18.5 arcmin 2 . Secondly, we obtained a 
dithered sequence of images of a central field of the galaxy 
providing deep photometry for a 4.5 x 4.5 arcmin 2 area. 
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Fig. 1. Map of the mosaic of 16 shallow fields observed 
in Fornax dSph, plotted on a DSS2-red image. The loca- 
tion of the deep field is shown as a grey box. 



We used an 8-points dithering pattern with shifts of up to 
10" from the central position. The observing parameters 
are given in Table[TJ which lists the number of images N of 
each sequence along with the Detector Integration Time 
(DIT) and the number of co-added integrations per image 
(NDIT). The on-target exposure times are 960s in J, H, K 
for the deep image, and 60s for the SHALLOW mosaic. The 
observation strategy is further illustrated in Fig. [TJ 

Given the non-negligible crowding of our Fornax field, 
we adopted an observing strategy based on offset sky im- 
ages alternated to on-target images. For the deep image, 
the "chopping" time interval varied between 120s for J, H 
and 60 s for the K band, while it was 60 s for all bands 
for the SHALLOW observations. The offset pattern was a 
dithered cross pattern for the deep image and a simple 
offset in declination for the SHALLOW scans. Although ex- 
pensive in terms of observing time, this strategy proved to 
give extremely good sky subtraction on scales comparable 
with the Point Spread Function (PSF) size, improving the 
quality of the photometry. 

Standard stars (including very red stars) from 
Persson et ail (jl998h were observed on a regular basis dur- 
ing the nights at airmasses comparable with those of target 
objects, to provide photometric calibration. For each stan- 
dard star, five images were taken, with the star located at 
the center of the detector and in the middle of the four 
quadrants of the frame. 
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Fig. 2. Residuals in the (a, 5) coordinates between the 
astrometric reference catalog (obtained from ESO WFI 
observations) and the SOFI catalog, after astrometric cal- 
ibration. The data refer to the whole wide-area mosaic, 
composed of 16 fields. The two histograms shows the Aa, 
Ad marginal distributions. The mean shifts are negligi- 
ble, 0'.'002 on both axes, while the standard deviation is 
O'.'IO and 0'.'12 on the right ascension and declination axis, 
respectively. 



2.2. Reduction and astrometry 

The redu ction steps were imple mented in as de- 



scribed in iMomanv et alj (|2003l ) . For the SHALLOW imag- 
ing, the basic observation and reduction unit was one 
"strip" consisting of 4 science and 4 sky frames. The sky 
images were scaled to a common median, after rejecting 
the highest and lowest pixels, and averaged to produce a 
master sky frame to be subtracted from all science images 
of the strip. The reduction steps were similar for the deep 
imaging, using the four offset sky frames closest in time to 
each science image. The sky-subtracted images were then 
corrected for bad pixels and flat-fielded. Illumination cor- 
rection frames (as well as bad pixels maps) available from 
the ESO Web pages were found suitable to this purpose. 
For the deep imaging, all individual images, after sky sub- 
traction, correction of bad pixels, and flat-fielding, were 
registered and averaged using the task imcombine with 
rejection of the brightest and faintest pixel for cleaning of 
cosmic ray hits. 



1 The Image Reduction and Analysis Facility (IRAF) 
software is provided by the National Optical Astronomy 
Observatories (NO AO), which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under 
contract to the National Science Foundation. 
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For the shallow scans, mosaicing is a complex task 
because of the modest overlap of the individual images. 
We resorted to absolute astrometry to effectively and ac- 
curately register the scans, match the photometry cat- 
alogs, and produce a single mosaic image. To this pur- 
pose, we used the IRAF package mscred (| Valded [l998h 
along with the wfpred script package designed at the 
Padua Observatory by two of us (LR, EVH) for reduc- 
tion of CCD mosaics. The first step in the mosaicing was 
the characterization of the astrometric properties of SOFI. 
This was done using a catalog of secondary astrometric 
standards in the central region of Fornax established from 
ESO-WFI observations (|Rizzi et al.ll2007h . This secondary 
catalog, in turn tie d to the USNO-A2.0 reference system 
( Monet et al.lll998h . is required to provide a sufficiently 
high surface density of stars needed to map distortions in 
the small field of SOFI. 

The field distortion of SOFI was actually found to be 
negligible. Once a distortion map was constructed, all the 
individual exposures were registered and resampled to a 
distortion-corrected coordinate grid with a common ref- 
erence World Coordinate System (WCS). As a result, a 
large-area mosaic image was reconstructed in each of the 
JHK filters. 

As a check of our astrometric calibration, we measured 
the right ascension and declination of the secondary stan- 
dard stars on the WCS-calibrated JHK frames. Figure [2] 
shows the difference between stellar coordinates on the 
SOFI wide-field catalog, and the a, S of the same stars in 
the reference optical catalog. This consistency test shows 
the internal precision of the astrometric calibration, which 
is the really important figure for image and catalog regis- 
tration. The absolute (systematic) accuracy of the coordi- 
nates is that of the reference catalog, estimated to be of 
the order 0.2 arcsec. The standard deviation of the resid- 
uals is not larger than C/12 on both coordinates. 



2.3. Stellar photometry and calibration 

Stellar photometry was obtained for the d eep and SHAL - 
low images using daophot ii/allstar ( Stetson! [l987l ). 
For the deep photometry, the PSFs were generated from 
a list of isolated stars on the coadded images. A Penny 
function with a quadratic dependence on position in the 
frame was adopted, to account for the elongation of the 
stellar profiles especially affecting the left side of the chip 
(this was traced to a misalignment of the Large Field ob- 
jective). The photometry was finally performed on the 
stacked images using ALLSTAR. For the SHALLOW images, 
simple aperture photometry with an aperture with radius 
equal to the FWHM of the stellar profile was found to 
yield the best photometric precision. For both catalogs, 
the pixel coordinates provided by daophot ii were con- 
verted to the equatorial system using the calculated WCS 
astrometric calibration. 
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Fig. 3. Comparison of the deep and SHALLOW photome- 
try of Fornax. 
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The photomet ric calibration techniq ues are similar to 
those adopted by iMomanv et al. (2003) and will be only 



Fig. 4. Comparison of shallow catalog with 2MASS 
photometry. Filled circles are data that passed a k<j re- 
jection and were used to compute the mean differences. 



briefly outlined here. Photometry of the standard stars 
through increasing apertures was obtained with the IRAF 
APPHOT task. For the "total" magnitude, we adopted a 
reference ape rture of 18 pixe l radiu s, close to the 10" aper- 
ture used bv IPersson et al. (1998). The 5 measurements 
of each standard star were averaged after checking their 
uniformity. The instrumental magnitudes were finally nor- 
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malized to 1 s exposure time and zero airmass; for a mag- 
nitude m ap , calculated for an observation with exposure 
time i cxp and airmass X , we defined the normalized mag- 



L cxp 

nitudc m as 



''ap 



2.5 log(texp) - K X X 



(1) 



where k\ are the mean atmospheric extinction coefficients: 
Kj = 0.08, kh = 0.03, and kk b = 0.05 were adopted from 
the ESO web page of SOFI. 

Given the small number of red stars observed each 
night, we used for calibration all data available for five ob- 
serving nights (including data previously obtained in Feb. 
2000 with the same SOFI setup). All the measurements 
were scaled to a common zero point and a first least square 
fit was done to compute the color terms of the calibration 
relations. Then, assuming fixed color terms, we measured 
the zero point variations through the run (in particular 
between the two nights of Nov. 10, 11) and found them 
to be very small, comparable with the measurements er- 
rors; a single zero point was therefore adopted for the run, 
with uncertainties 0.02 mag in J, H , and K. The resulting 
calibration relations are: 



J-j = -0.016 (J- H) + 23.118 
H-h = +0.001 (J — H) + 22.902 
H — h = -0.002 (H-K)+ 22.902 
K - k s = +0.021 (J - K) + 22.337 



(2) 



which are consistent with the calibration presented in 
the ESO web page . Note that our magnitudes on the 
Persson et al.l (|l998T ) K system are expected t o show vir- 
tually no offset relative to the 2MASS system ([Carpenter 
2001 ). The c o lor te rms measured with respect to the 



Persson et al.l (J1998J) JHK standard stars are negligible. 

The photometric catalogs of Fornax stars were cali- 
brated using the Eqs. ([2]) after magnitude-scale and aper- 
ture correction. This was based on APPHOT large-aperture 
photometry and growth-curve analysis of a few relatively 
bright, isolated stars in the deep and SHALLOW fields. The 
uncertainties on aperture correction are of the order 0.03 
mag, yielding total calibration uncertainties 0.04 mag. 

In order to test the photometric calibration and rule 
out the presence of photometric bias in the shallow pho- 
tometry, we compare the zero points of the deep and 
shallow photometry in Fig. [31 For the stars in common 
the photometry shows an excellent internal consistency 
down to K « 18, with zero point differences of the or- 
der 0.01 mag. Figure [4] shows a compa rison of our SHAL- 
low catalog with 2MASS photometry (jCutri et al.ll2003h 
on the same area. We have selected only 2MASS mea- 
surements with S/N > 10. The systematic differences are 
negligible considering the errors in the 2MASS catalog and 
the aperture correction uncertainties. 

The photometry lists obtained for the individual point- 
ings were merged into a single shallow catalog. The stars 
were matched on the basis of right ascension and declina- 
tion, using an error box of 1'.'5. Finally, the coordinates of 




mag 



Fig. 5. Completeness and photometric errors from artifi- 
cial star experiments for the SHALLOW photometry. Upper 
panel: Completeness calculated as the number of artificial 
stars recovered as a function of observed magnitude. Lower 
panel: photometric errors, i.e. the standard deviation of 
the differences between input and output magnitudes. 

stars in the two catalogs were matched to those of BVI 
photometry (select ing objects with V < 23) from ESO 
WFI observations (jRizzi et al.ll2007j ). thus producing a fi- 
nal list with BVI J UK magnitudes. The SHALLOW and 
deep catalogs of near-infrared photometry of Fornax stars 
are presented in Tables [2] and [3] (the entire catalogs are 
made available electronically at the CDS). 



2.4. Artificial star experiments 

The completeness of our photometric catalogues was eval- 
uated from artificial star experiments. We performed 20 
test runs by adding 400 stars to the scientific frames in 
each run. Since the crowding is uniformly low over the 16 
SHALLOW pointings, with only about 1000 stars per point- 
ing on average, the experiments were limited to one of 
the shallow pointings. The input magnitudes and col- 
ors of the artificial stars were chosen along a sequence 
corresponding to the Fornax RGB. The stars were placed 
on a grid of equilateral triangles with a side of 40 pixels 
(~ ll'/5), much larger than the stellar PSF. In each ex- 
periment the grid was then randomly shifted in order to 
uniformly cover all the frame. The results from artificial 
star experiments on the SHALLOW and deep photometry 
are shown in Fig. [5] and Fig. [6J respectively. We note that 
the deep data are comparatively rather noisy in the H 
band with respect to the J and K bands. 

Another well-known observational effect is the fact 
that the mean difference between input and output mag- 
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Table 2. The SHALLOW near-infrared catalog of Fornax stars over a 18.'5 x 18'5 area. A few lines are shown here for 
guidance regarding its form and content, while the full catalog is available from the CDS. 
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2 40 22.48 
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Table 3. The first rows of the deep near-infrared catalog of stars in a 4f5 x 4.'5 central region of Fornax. The entire 
catalog is available only electronically at the CDS. 
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Fig. 6. The results of artificial star experiments for the 
deep photometry. 

nitudes is g enerally biased tow ards brighter output magni- 
tudes (e.g., Gallart et al .II 1 9 96h . In our experiments for the 
shallow data, the mean difference of input and output 
magnitudes resulted to be less than 0.05 mag for K < 18.2. 
This is also an upper limit for the (J—K) color shift down 
to this limiting magnitude. However, the color bias is neg- 
ligible (< 0.01 mag) in the first 2 mag below the RGB tip 
used for the metallicity measurements. 



3. Color-magnitude diagrams 

The near-infrared and optical-infrared color-magnitude di- 
agrams (CMDs) of Fornax dSph in the SHALLOW sample 



are presented in Fig. [JJ Given the better image qual- 
ity and spatial resolution of our optical photometry, we 
used the ALLSTAR shape parameter sharp for the V mea- 
surements to remove bad and non-stellar objects from our 
photometric catalogues. Only objects with —0.3 < sharp 
< 0.3 were selected. 

The most noteworthy feature in the CMD in Fig.[7Jis 
the well defined sequence of intermediate-age AGB stars. 
A fraction of AGB stars have red colors typical of car- 
bon (C) stars, and a few stars show extremely red col- 
ors possibly indicating dust-enshrouded AGB stars. The 
properties of near-infrared selected AGB stars in Fornax 
will be discussed in more detail in a separate paper. We 
only remark here the dramatic change in their distribu- 
tion when optical and near-infrared bands are used. Note, 
for example, that the redder AGB stars become progres- 
sively fainter in the bluer optical bands, so that even the 
brightest (in terms of bolomctric luminosity) arc missed in 
optical CMDs. For the reddest AGB stars the (B-V) col- 
ors saturate. Thus, selection in the near-infrared appears 
to be extremely important to investigate the evolution of 
luminous AGB stars. Also, an advantage of optical-near- 
infrared colors is the improved discrimination against field 
contamination by foreground stars and background galax- 
ies, as is evident, for example, when comparing the V, 
(B - V) and K, (J - K) diagrams in Fig. 

The shallow, large-area catalog provides the statistics 
for studying the AGB stars, the RGB tip, and the stel- 
lar populations of Fornax. The fainter limiting magnitude 
and better precision of our deep photometry of a smaller 
central region allows us to characterize the RGB down to 
the so-called "AGB bump" and the red clump. Our deep 
color-magnitude diagram is presented in Fig. together 
with the RGB fiducial li nes of Milky Way globular clus- 
ters (IValenti et al. 2004a). Thi s diagr am is similar to that 
obtained bv IPietrzvhski et all (|2003l ) with ISAAC at the 
VLT. Note that the RGB of Fornax is relatively thin in 
this small central field. It lies between the fiducial RGB 
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Fig. 7. The color-magnitude diagrams of Fornax from the SHALLOW catalog, obtained by combining optical and 
near-infrared photometry. Stars on the extended AGB, chosen redder than (J — K) = 1 and brighter t han K = 14.5, 
are sh own as filled circles. Su perimposed on the K, (J — K) diagram are theoretical isochrones from iGirardi et alJ 



( 20021 : G02, dashed lines) and Pietrinferni et al. ( 2004 : P04, solid lines). The isochrone age is 7 Gyr and the metal 
abundances are the closest to the metallicity of Fornax for each model set: Z = 0.001 and 0.004 for G02, Z = 0.001 
and 0.002 for P04. Typical photometric errors as derived from artificial star experiments are also shown as crosses in 
the upper-left panel. 



lines of M 4 ([Fe/H] = -1.19 on the ICarretta fc Gratton 
19971 abundance scale) and M 107 ([Fe/H] = -0.87). This 
suggests a metallicity [Fe/H] ~ -1 on the same scale. 

4. Luminosity function and distance 

Using the deep data shown in Fig. El we constructed a 
luminosity function (LF) extended down to the red clump 
of Fornax (Fig. [9]). The LF was corrected by taking into 
account the incompleteness factor derived in Sect. 12.41 
Our sample of RGB and He-burning stars was selected 



as shown in the inset. Using this LF, new estimates of the 
distance to Fornax have been obtained based on the mean 
A'-band luminosity of the red clump and the K magnitude 
of the RGB tip. 



4.1. The red clump 

Several authors have pointed out that distance measure- 
ments based on the magnitude of the RC are most re- 
liable in the near-infrared, because of the smaller dc- 
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47Tuc -0.70 
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Fig. 8. The near-infrared CMD of Fornax derived from 
our deep photometry, showing the RGB and the promi- 
nent "red clump" of intermediate-age He-burning stars. 
Superposed on the R GB are the fiducial lin es of Galactic 
globular clusters from Valenti et al. (l2004al). with metal- 
licities on the scale of Carretta fc Grattonl ( 199?t) . Typical 
photometric errors as derived from artificial star experi- 
ments are shown as crosses. 



pendence of the RC mean luminosity on mctallicity and 
age in the A"-band, along with the reduce d reddenin, 

Alvesl I200 



relative to the optical wavelengths (e.g 



Grocholski fc Saraiedinil l2002t IPietrzvhski et all l2003h 
Pietrzvhski et all (|2003l) employed VLT-ISAAC photome 
try of RC stars to measure the Fornax distance. An inde- 
pendent distance estimate based on our own measurement 
of Arc is presented here, and the two determinations com- 
pared. 

In order to reduce the dependence from the bin choice, 
the LF was constructed as a multibin histogram by av- 
eraging 10 LFs with a fixed 0.1 mag bin and starting 
points shifted in steps of 0.01 mag. A mean level Arc = 
19.20 ± 0.02 was measured following the standard proce- 
dure, i.e. by fitting the sum of a Gaussian and a polyno- 
mial to the magnitude distribution of stars in the color 
range 0.30 < J — K < 0.85 (see, e.g., Pietrzvhski et al.l 
20031) . The mean Arc value and its standard deviation 
are the results of 5000 experiments with bootstrap resam- 
pling of the luminosity function. Based on our artificial 
star experiments, any magnitude shift due to a photomet- 
ric bias is negligible (~ 0.01 mag) for the deep data at 
the RC level. Other sources of error are considered in the 
following, in addition to the 0.02 mag fitting error, to eval- 
uate the total uncertainty in the distance to Fornax dSph. 
In order to compute th e distance, our ma gnitudes 



Fig. 9. The A-band luminosity function of evolved stars 
in Fornax, corrected for completeness effects. This LF, de- 
rived from the DEEP catalog, shows the location of the red 
clump and the AGB bump. The inset illustrates our sam- 
ple selection. A histogram with a 0.1 mag bin is shown 
as a dotted line, while the solid line represents a multibin 
histogram. 



be transformed onto the pho tometr i c syst em used in the 
RC luminosity calibration of lAlved (l2000l) . Since our A- 
band photometry agrees very well with the 2MASS system 
(Sect. l2~3l see also lCarpenterl2001 ) , we simply adop t ed the 
transformation from 2MASS to the lBessell fc Brettl (|l988l ) 
syste m, Abb = Admas s + 0.044. The l atter is very close 
to the Koornneej (Il983l) system used by Alvesl ( 2000l) (see 



Carpenterll200ll IGrocholski fc Saraiedinill2002l ) 
We therefore applied the relation 



(to - M) = (Arc + 0.044) - M 



A" 



.4 



K 



AM 



K 



(3) 



where the A-band luminosity of the red clump in 
the so lar ne i ghbor hood is Mk = —1.61 according 
to the Alvesl (2000h calibrat ion (a value confirmed by 
IGrocholski fc Saraiedini 2002 from a sample of 14 open 
clusters), and AMk is a population correction term. This 
correction accounts for the different stellar cont ent of 
Forna x and the local Galactic RC on which the Alvesl 
(2000) calibration is based. The population correct i on wa s 
calculated using the precepts of Salaris fc Girardi ( 20021 ). 



(tied to the LCO system of iPersson et al.lll998l) need to 



the age-metallicity relati on oflPont et al.l (120041) . and the 
star-formation history bv lTolstov et all (|2003l) . and found 
to be AMk = —0.10 (the RC in Fornax being fainter) 
Adopting the extinction law of iRieke fc Lebofskvl (Il985h 
and E(B — V) = 0.03, the resulting distance is (m — M)o = 
20.84 + AM K =20.74 ± 0.11, where the uncertainty in- 
cludes the statistical error on the RC location, the pho- 
tometric zero-point error, and a photometric error of 0.1 
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mag a t the level of the RC (see Fig. [5]). iPietrzviiski et al 
( 20031) . considering the population correction negligible, 
found a distance modulus (to — M) — 20.86, which dif- 
fers from our value only for the correction term, the -/^rc 
determination being in perfect agreement. 

4.2. The AGB bump 

In the LF presented in Fig. [9] another bump is clearly seen. 
Its K band magnitude was measured by fitting also in this 
case the sum of a polynomial and a Gaussian function 
and found to be K = 17.99 ± 0.03 with a bootstrap tech- 
nique similar to that used for the red clump. This feature 
is identified with the AGB bump, which is the signature 
in the CMD of the beginning of the AGB phase. At the 
beginning of this evolutionary phase the increase of lumi- 
nosity of a star is slower than in subseq uent AGB phase, 
and t hus a bump in the LF is produced dCastellani et al 



19911 : lGallartlll998l : lAlves fc Saraiedinilll999D 

We also explored the possible identific ation of this 
feature with the RGB bump, using the IValenti et al 



1 2004bJ ) calibration of the RGB bump. Assuming a metal- 
licity [M/H1= —1-6 for the old popul ation of Fornax 
( Bersier fc Woodll2002l ; [Greco et all2006f) . we derived that 
the expec ted RGB bump magn itude is K = 18.52. As dis- 
cussed bv lSaviane et al .1 ( 2000a ). the RGB bump is fainter 
that the bump in Fig. [5] and not visible because it is too 
close to the overwhelming red clump. We note that the 
metallicity assumed here, [M/H]= —1.6, is appropriate 
for the old population of Fornax and a lower limit to the 
actual mean metallicity of intermediate age stellar pop- 
ulations (see fol lowing section s ). The RGB bump magni- 
tude derived bv lValenti et al. ( 2004bl ) becomes fainter at 
increasing metallicity, and therefore our estimate of the 
RGB bump magnitude is to be considered a lower limit. 
We therefore identify the feature observed at K = 17.99 
with the AGB bump. 

4.3. The RGB tip 

Figure [10] shows a close-up view of the brighter part of the 
RGB luminosity function in Fornax, obtained by select- 
ing red giant stars from the wide field SHALLOW catalog 
This catalog was chosen because of the better statistics, 
and, indeed, the RGB cutoff appears very well defined. 
We obtained an objective estimate of the magnitude of 
the RGB tip by fitting the LF with the convolution of 
a step function with a Gauss ian kernel representat ive of 
the measurement errors, as in lMomanv et all ( 2002h . The 
function is composed of a constant value brighter than the 
RGB cutoff, and a power law (on a logTV scale) below the 
RGB cutoff. Since the tip is located at K ~ 14.6 where 
our photometry is complete (see Fig. [5]) , no completeness 
correction was applied to the LF. As before, a multibin 
LF was used, obtained by averaging 10 LFs with a fixed 
0.1 mag bin and intervals shifted by 0.01 mag. Using this 
procedure the RGB tip was detected at K = 14.59 ± 0.03. 




0.5 



[M/H] 




(V-K) 

Fig. 10. The luminosity function of the upper RGB, from 
the shallow catalog. The inset shows how RGB stars 
have been selected. The RGB cutoff was measured by fit- 
ting an error-convolved step function (green long-dashed 
line) to a multibin histogram (continuous line). The dot- 
ted line shows an example of a simple LF with 0.1 mag 
bins. 



The uncertainty includes the fitting error and the error 
associated with the binning of the LF. The internal mea- 
surement error at the level of the RGB tip gives a minor 
contribution (~ 0.01 mag). 

An independent measurement of the RGB tip was also 
obtained u sing the Maximum L ikelihood Algorithm de- 
scribed bv iMakarov et all (|2006l ). The TRGB was found 
at K = 14.63 ± 0.03. Since the two independent measure- 



ments agree within the errors, we adopt the mean value 
K = 14.61 ± 0.02 (random) ±0.03 (systematic) as our fi- 
nal measure of the tip, where the systematic uncertainty 
reflects the photometric zero point error (Sect. H]). 

The age and metallicity dependence of the RGB tip 
is larger in K than in the /-band, so that the applica- 
tion of th e RGB tip method is more uncertain in this 
case (e.g., ISalaris fe Girardill2005l) . Intermediate-age stars 
show a fainter RGB tip than old stars, with a difference 
that can be as large as ~ 0.1 mag in K in the age inter- 
val 4-13 Gyr. On the other hand, if younger populations 
become more metal-rich as a result of galactic chemical 
evolution, their RGB tip becomes brighter, because the 
TRGB K luminosity rises with increasing metallicity. 

Since the bulk of stellar populations in Fornax is 
of intermediate age, a population correction is required 
to compute the distance. To this aim, we modeled the 
combined effects of age and metallicity on the RGB 
tip by constructing a synthetic CMD containing 100 000 
stars in the upper RGB. Our simulations arc based 
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on the ZV AR code (IBertelli et al.l Il992f ), the Padova 
isochrones ( Girardi et alj 2002 ), and adopt the chemical 



evolution history by iPont et all ( 2004 ) and the SFH from 
Tolstoy et ail (|2003h . We obtained the population correc- 
tion for the luminosity of the RGB tip by measuring the 
cutoff in the simulated CMD for (i) the full stellar popula- 
tion mix (representative of the Fornax RGB), and (ii) only 
stars older that 10 Gyr. As a consequence of the adopted 
metal enrichment history, the TRGB in the overall syn- 
thetic CMD turns out to be brighter than that for old 
stars alone. The difference in the TRGB K magnitude is 
AM K = M$ d - Mf = 0.22 mag. 

By taking this population correction into account, 
the distance modul us of Fornax was calculated using 
Valenti et all (|2004bl ) empirical calibration of the M£ RGB 
as a function of metallicity, based on Galactic globular 
clusters: 



M 



TRGB 



-6.92 - 0.62 [M/H] 



(4) 



where [M/H] is an estimate of the average metal abun- 
dance. The r.m.s. uncertainty of this calibration is 0.16 
mag. We then apply Eq. [4] to the old stellar population in 
Fornax: 



(m - M) = (K 



TRGB 



AM K ) - A K - Mk 



TRGB 



(5) 



where K TKGB i s the measured TRGB level and A K 
is the extinction. Assuming a low metallicity [M/H] = 
— 1.60 ±0.15 appropria t e for the old popul ation of Fornax 
( Bersier k Wood! 120021 : fGreco et al.ll200i ). the corrected 
distance modulus is (m — M)o = 20.75 ± 0.19, where the 



error is dominated by the uncertainty on the I Valenti et al 
( 200411 calibration. 

It is interesting that, once the stellar content is prop- 
erly taken into account, the if-band RGB tip provides 
an estimate of the distance to Fornax in good agreement 
with our determination from the RC and other estimates 
in the literature. However, the uncertainty remains high, 
given the strong dependence of the TRGB on age and 
metallicity and the error on the calibration relation. 

In summary, the distance modulus derived here from 
near-infrared data and methods ap pears to confirm those 
measured by Saviane et al.l ( 2000al ) from the I magnitude 
of the RGB tip, (to — M) Q = 20.70, and from the mean 
magnitude of old HB stars, (m — M) = 20.76. New dis- 
tance measurements b ased on wid e -field optical observa- 
tions are presented by iRizzi et al.l ( 2007 ). where the dif- 
ferent results in the literature are compared and discussed 
in some detail. 

5. Metallicity 

5.1. Mean metallicity 

The mean metallicity of the stellar populations that 
make up the Fornax RGB were estimated by compar- 
ing the near-infrared and optical-infrared colors of red gi- 
ants to the RGB fiducial lines of Galactic globu l ar clus- 
ters of known metal abundance. IValenti et al.l (|2004a[ ) 



Table 4. The mean metallicity of Fornax from 

photometric indice s of IValenti et alJ (j2004al ) and 

Kuchinski k Frogell (| 19951 ). 



index 




[Fe/H] gc 


[M/H] 


(J~K)o K= 


-5.5 


-0.98 ±0.18 


-0.81 ±0.17 


(J - K)^= 


-5.0 


-0.97 ±0.20 


-0.83 ±0.19 


(J - K)™«= 


-4.0 


-0.98 ±0.19 


-0.86 ±0.18 


(J - K)?*= 


-3.0 


-1.01 ±0.15 


-0.79 ±0.14 


(V - K)™«= 


-5.5 


-1.13 ±0.15 


-0.96 ±0.16 


(V - K)^ 


-5.0 


-1.23 ±0.27 


-1.02 ±0.25 


(V - K)^= 


-4.0 


-1.37 ±0.39 


-0.94 ±0.35 


(V - K)™«= 


-3.0 


-1.28 ±0.38 


-1.04 ±0.39 


VFO slopeaGB 


-1.20 


-1.03 


KF slope R ,GB 




-1.16 





calibrated the RGB colors at fixed ii'-band luminosi- 
ties in the near-infrared CMDs of Milky Way globular 
clusters as a function of metallicity. They give color- 
metallicity relations for (V — K)o and (J — K)q at Mk = 
-5.5, - 5.0, -4.0, -3.0 against the [Fe/H] of GGCs on the 
scale of ICarretta k Grattonl (|l997l ). 

They also provide calibrations against [M/H], a mean 
metallicity measuring the abundance of all heavy ele- 
ments. This parameter is particularly important to es- 
timate the metallicities of dwarf spheroidal galaxies by 
comparison with the photometric properties of Milky 
Way globular clusters. These objects are known to show 
non-solar abundance patterns, with an overabundance 

of a-clcmcnts relative to iron th at is a function of 

the cluster metallicity ( Pritzl et al. 20051 and references 
therein). In contrast, dwarf spheroidal g alaxies tend to 



have [ a/Fel ratios c loser to solar (e.g., IShetrone et al 



2003 ). Salaris et al. ( 19931 ) have shown that the color of 



red giant stars is driven by the overall metal abundance 
rather than the Fe abundance. Thus, the iron [Fe/H] scale 
of Galactic globular clusters is not immediately applicable 
to dwarf galaxies. Instead, the [M/H] provides a suitable 
parameter for comparing stellar systems with different 
abundance patterns and rank them againsts Milky Way 
globular clusters. For the sake of comparison with previ- 
ous works, we provide here both the [Fe/H] and [M/H] 
rankings, but recommend the [M/H] as the most appro- 
priate values. 

The mean [Fc/H] and [M/H] values of Forn ax RGB 
stars computed using the IValenti et al.l (|2004al ) calibra- 
tions arc presented in Table [4j together with errors 
calculated from color uncertainties by error propaga- 
tion. The average values are [Fe/H] y _ K = —1.22 ± 
0.12 and [Fe/H] r „ = -0.99 ± 0.09 (on the scale of 



ICarretta k Grattonlll997l ) . Using the calibrations against 
[M/H], wc obtain a mean metallicity [M/H] y Ar = 
-0.98±0.12and [M/R]j_ K = -0.82±0.08 Alternatively, 
a robust m etallicity indicator is repres ented by the RGB 



slope. The Kuchinski k Frogel (11995) ca li bratio n yields 
[Fe/H] zw = -1.28 on the IZinn k West! (|l984l ZW84) 
scale, corresponding to [Fc/H] : 



-1.16 on the CG97 
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Table 5. Galactic globular clusters used a s RGB metal- 
licity temp lates. We list metallic ities on the Zinn fc Westl 
( 1984) and lCarretta et all (|200ll ) [Fe/H] scales, as well as 
"global" metallicities [M/H] (see text). 



cluster 


[Fc/H] zw 


[Fe/H] CG 


[M/H] 


M15 


-2.17 


-2.12 


-1.91 


M68 


-2.09 


-1.99 


-1.81 


M55 


-1.82 


-1.61 


-1.41 


M13 


-1.65 


-1.39 


-1.18 


NGC 288 


-1.40 


-1.07 


-0.85 


NGC 362 


-1.27 


-1.15 


-0.99 


47Tuc 


-0.71 


-0.70 


-0.59 


NGC 6624 


-0.63 


-0.63 


-0.48 


NGC 6553 


-0.29 


-0.06 


-0.05 


NGC 6528 


-0.33 


+0.07 


+0.05 



scale (using the conversion of ICarretta et al.ll200lh . This 
val ue is confirmed by the recent re-calibration of the slope 
by IValenti et all (|2004al ). yielding [Fe/H] = -1.20 or 
[M/H] = —1.03, in good agreement with other methods. 

5.2. Metallicity distribution 

The metallicity of RGB stars can also be estimated by 
interpolating the colors of individual stars across a grid of 
empirical RGB templates. In this way, we also obtain an 
observable related to the metallicity distribution function 
(MDF) of the red gia nt stars. The method used here (see 
ISaviane et al.l 2000b ) consists in building up a family of 
hyperbolae that best fit the RGB fiducial lines of GGCs 
of known metallicity, and then use interpolation over that 
family of curves to derive the metallicity of RGB stars. A 
previous application of the method t o the [ V — K), Mr 



plane can be found in IZoccali et alj (|2003l ) , where it is 



discussed in some detail for stars in the Galactic bulge. 

Table \E\ lists the reference globular clusters with 
their adopted metallicities. In the present implement a- 
tion of the method w e have employed [M/H] values from 
Valenti et al.l (12004a ). with the updated [Fe/H] CG scale of 
Carretta et al.l (|200ll ). 

The resulting "photometric MDF" is shown in Fig. [TT1 
along with an illustration of the interpolation method. 
We also estimated the statistical uncertainty of our pho- 
tometric metallicity determinations using a Monte Carlo 
approach. A synthetic CMD was generated by randomly 
choosing 10 000 stars along the best fit line in Fig. [11] 
Errors were added to the synthetic magnitudes accord- 
ing to a Gaussian distribution with standard deviations 
0.01 in V and according to the results of artificial star 
experiments in K (see Sect. 12. 4p . We then used the same 
algorithm applied to the Fornax CMD to retrieve individ- 
ual metallicities. The recovered metallicities have a nearly 
Gaussian distribution with a dispersion a = 0.06 dcx. This 
scatter was taken as the statistical uncertainty associated 
to our individual metallicities. We note that the photo- 
metric completeness in the region of the CMD used to 




0.5 



[M/H] 




(V-K) 

Fig. 11. Derivation of a photometric Metallicity 
Distribution Function of RGB stars in Fornax using 
(V — K) colors. Upper panel: the photometric metallicity 
distribution function of red giants in Fornax. The solid 
line is a fit of the sum of two Gaussian functions [dashed 
lines) to the data. The method is illustrated in the 
lower panel: the stars on the upper RGB of Fornax are 
interpolated across analytical fits [dotted lines) to the 
RGB fiducial lines of template Galactic globular clusters. 
The yellow line marks the best-fit mean metallicity. 

compute the MDF is 100%, so that no completeness cor- 
rection is needed. 

The optical-near-IR data confirm the extended 
metallicity di s tributi on of Fornax stars suggested by 
Saviane et al.l (|2000al K The distribution is formally mod- 
eled by the sum of two Gaussians, with a main peak at 
[M/H] ~ -1.02 and a secondary peak at [M/H] ~ -1.28. 
The mode corresponds to about [Fe/H] ~ — 1.2 and —1.3 
on the CG97 and ZW84 scales, respectively. The standard 
deviation of the Gaussian corresponding to the main peak 
is 0.11 dex, which is much larger than the scatter of 0.06 
dex due to photometric error. We remind that our "metal- 
licity distribution" is representative of the MDF only for 
old stellar populations with age comparable to the old 
globular clusters in the Milky Way. This is certainly not 
true for Fornax stars , whose mean age is about 5-6 Gyr 



(jSaviane et al.ll2000al ). 

In all cases, however, this distribution represents an 
important observable that models of galactic evolution 
should be able to reproduce. 

5.3. Age correction 

Since the bulk of the stellar populations in Fornax is 
younger than stars in GGCs, the red giants stars in Fornax 
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Fig. 12. Contours of constant (V — K) color at Mk = 
—5.0 as a function of age and [M/H] (dotted line s), us- 
ing the model isochrones of IPietrinferni et all (|2004T ). This 
plot allows a quantification of the age-metallicity degener- 
acy in the mctallicity estimates based on optical-infrared 
colors. The red solid line corresponds to the color of the 
RGB of Fornax. An arrow shows the change in metallicity 
at this fixed color between 12.5 Gyr old globular clusters 
and a 7.5 Gyr old galaxy (representative of Fornax dSph). 
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Fig. 13. Comparison of metallicities of individual stars 
inferred from near-infrared colors and metal abundances 
derived from the strengths of Call triple t lines . Lower 
panel: spectrosc opic data from 



Pont et al 



panel: data from IBattaglia et al 



(|2004D ; upper 
(|2006h . The solid line in- 
dicates one-to-one correlation, while a dashed line shows 
the shift expected by taking into account the mean age of 
stellar populations in Fornax. Error crosses represent the 
mean uncertainties of spectroscopic measures, and the 3cr 
error of our photometric determination. 



are on average slightly bluer than globular cluster stars of 
the same metallicity. As a consequence, the metallicity 
obtained from the mean RGB color is systematically un- 
derestimated. Using optical-near-infrared colors, however, 
this "age-metallicity degeneracy" is much reduced with re- 
spect to optical colors. This is illustrated in Fig.[T21 where 
we show the effects of age and metalli city variations on 



RGB colors at Mk = —5.0 using the IPietrinferni et al 



( 2004J) isochrones. In this figure, contour lines of equal 
color are nearly vertical. This means that the effects of 
a change in age (or those of an age spread) on the color 
shift and the color dispersion are much smaller than those 
produced by mctallicity variations. 

We have used this plot to estimate a mean correction to 
metallicity, by assuming that the age of GGCs and Fornax 
dSph arc 12.5 and 7.5 Gyr, respectively. The correction 
is the difference in mctallicity needed to keep the color 
constant while moving from 12.5 Gyr to 7.5 Gyrs. We 
found a correction A [M/H] ~ 0.15, yielding for Fornax 
an age-corrected metallicity [M/H] = —0.87 in terms of 
global metallicity (we consider the mode of the photomet- 
ric MDF). 

This result is i n good agreement wit h previous results 
in the literature. ISaviane et alj ( 2000a ) derived a mean 
age-corrected metallicity [Fe/H] = —1.0 from the (V — I) 



color of the RGB, on the ZW84 scale (corresponding to 
[Fc/H] CG « -0.9). 

5.4. Comparison with spectroscopy 

Our mean metallicity agrees well with the results of low- 
resolution spectroscopic analyses of RGB stars, yielding 
a mean metallicity [Fe/H] = -1.0 and [Fe/H] = -0.9 
(|Tolstov et al.ll200ll ; IPont et al.|[200l on the CG97 scale. 
This value has been conf irmed by high-resol ution spec- 
troscopy of (two) stars bv lTolstov et al. ( 2003 ). Since the 
extended star for mation of Fornax leads to an abundance 
ratio [a/Fe] ~ ( Shetrone et al. 2003T ). we can assume 
that our [M/H] values are directly comparable with the 
spectroscopic measurements. 

The relatively large overlap between our optical-near- 
IR photometry and the spectroscopic samples of Fornax 
RGB star s with metallicities derived by Ca H tri plet spec- 
troscopy ( Pont et ail2004 ; Battaglia et al.|[2006h allows us 
a direct comparison of metal abundances on a star-by-star 
basis. This comparison is especially interesting to assess 
the reliability of photometric MDF determinations for all 
systems that are too distant for (even low-resolution) spec- 
troscopy. 
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the mode of the metallicity distribution with the results 
of spectroscopy. 

Clearly, the first explanation for the discrepancy at the 
metal-rich end is the age-metallicity degeneracy in RGB 
star colors that even optical-infrared color indices can- 
not completely overcome. Assuming that the Can triplet 
methods provides the correct metallicities, the behavior 
noticed in Fig. [13] and [14] appears to be consistent with 
a late metal enrichment scenario suggested by the cited 
spectroscopic studies. Indeed, the redder optical-infrared 
colors of a young population of metal-rich red giants are 
compensated for by the younger age. However, spectro- 
scopic abundances from Can triplet line strengths are 
also somewhat uncertain, especially in the high-mctallicity 
regime (where an extrapolation may be needed) and when 
the Call triplet calibration, which is based on globular 
cluster stars, is applied to the spectra of young (1-4 Gyr 
old) stars. Accurate abundance measurements from high- 
resolution spectroscopy may be useful to definitively clar- 
ify the issue. 



Fig. 14. Comparison of the metallicity distribution in- 
ferred from the optical-infrared colors of RGB stars (this 
paper, upper panel) a nd the MDF obt a ined from Can 
triple t spectroscopy ( Pont et al. 2004 ; Battaglia et all 
20061 lower panel) . The shaded part of the histogram from 



Pont et al 



2004h ([Fe/H] > -0.8) contains extrapolated 
metallicities. A mean age correction, resulting in a shift of 
0.15 dex, has been applied to our histogram. 

Figure [T3l plots the metallicities of individual stars de- 
rived from (V—K) colors against those estimat ed from the 



equiva lent widths of the Call triplet lines by | Ppnt et al . 
(|2004f) (lower panel) and iBattaglia et~ai1 (120061) (upper 



panel). For the comparison with Pont et al. ( 20041 ) results 
we used only spectra with noise below a given thresh- 
old (F. Pont, priv. comm.). The spectroscopic values for 
metal-rich stars are those corrected by the authors by 
comparison with high-metallicity stars in the LMC (see 
Pont et al.ll2004 their section 3). While a few stars ap- 
pear too red or too blue for their spectroscopic metallic- 
ity, the general trend is that of an overall correlation. A 
discrepancy is apparent at the metal-rich end, where the 
metallicity estimates from photometry appear to saturate. 
A similar trend, with a worse cor r elatio n, is noticed in the 
comparison with Battaglia et al. ( 20061 ) data, where there 
is a large excess of metal-rich stars with respect to the 
photometric estimates. 

The metallicity distribution of Fornax RGB stars is 
shown in Fig. Q3J together with the MDFs obtained from 
Call spectroscopy. In this case, the distribution of photo- 
metric [M/H] values in Fig. [Ml was corrected by 0.15 dex 
toward higher metallicities to take into account the fact 
that intermediate age stars (the bulk of Fornax stars) are 
bluer than the GGC template stars of the same metallic- 
ity. Although this correction clearly represents a first order 
approximation, it is interesting to note the agreement in 



6. Summary and conclusions 

We have presented a near-infrared photometry of the stars 
in the Fornax dwarf spheroidal galaxy. Our study provides 
color-magnitude diagrams and photometric catalogs of red 
giant and AGB stars in Fornax over a 18f 5 x 18'5 area, 
and deep photometry over a 4.'5 x 4.'5 central field. The 
main results are the following: 

— From stars on the red giant branch and the RC, we 
have obtained independent estimates of the distance to 
Fornax based on mean K magnitude of the red clump 
and the RGB tip, which take into account the mean age 
of the stellar populations in Fornax. The average value 
obtained from the two methods is (to — M)o=20.74± 
0.11, in excellent agreement wit h previous authors, and 
in particular with the results of Saviane et alJ ( 2000a[ ) 
from optical photometry. 

— The (V — K) color distribution of RGB stars has been 
used to infer the mean metallicity and metallicity dis- 
tribution of red giant stars in Fornax, taking advantage 
of the reduced dependence of (V — K) colors of RGB 
stars from the age of the stellar population. The av- 
erage metallicity was found to be [M/H] ~ —0.9. This 
compares well with the values recently ob t ained from 



spectr oscop y by Tolstoy et a l . d200ll . l2003f) , IPont et al 
(|2004f) . and iBattaglia et al ' " 



( 20061) . 



The metallicity distribution is consistent with that 
obtained from spectroscopy up to the metallicity of 
47Tuc ([Fc/H] ~ —0.7). However, there is a clear 
discrepancy between the MDFs derived from near- 
infrared co lors and sp e ctrosc opy near the metal-rich 
end, where IPont et all (|2004h found of a tail of stars 
with metallicity up to almost solar. This discrepancy 
could be caused by the effects of the age-metallicity de- 
generacy, which cannot be completely corrected even 
using ncar-IR photometry, but the extrapolation used 
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bv lPont et al. ( 2004 ) to derive their metallicities from 
the measurements of Call lines, as noted by the au- 
thors could also have some effect. More observations 
are needed to solve this discrepancy and establish the 
upper end of metal enrichment in Fornax. 
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